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G^ ' We study the constraints on the CP-violating soft-breaking phases in the minimal supersym- 

metric standard model using the limits on the chromoelectric dipole moment of the strange 
Y^. quark extracted from the neutron EDM experiment. Our investigation shows that the phase 

Qh! mediated by the gluino exchange diagram has to be very small, (f) < 8 ■ 10^^, for the com- 

(-H I mon supersymmetric mass of the order of 100 GeV. Then, solving the renormalization group 

equations analytically by iterations, we calculate the electric dipole moment of the neutron 
in the MSSM with CP-conserving soft-breaking parameters for the case of three and four 
generations. For the three-generation case we resolve the apparent discrepancies between 
order-of-magnitude estimates and numerical calculations existing in the literature. In this 
case the EDM of the neutron does not exceed 10~^^e ■ cm. For the four-generation case we 
show that there is a significant enhancement which renders the EDM of the neutron at a 
measurable level of 10~^^e • cm. 
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1 Introduction 

The current experimental limit on the electric dipole moment (EDM) of the neutron |[l[], 
(In /e < 10^^^ cm, exceeds the realistic Standard Model prediction for this quantity by seven 
orders of magnitude. This gap between theory and experiment provides an excellent oppor- 
tunity to limit a new CP- violating physics beyond the Standard Model (SM). The purpose 
of this work is to consider the electric dipole moment of the neutron in the Minimal Super- 
symmetric Standard Model (MSSM) with three and four generations of quarks. 

The Minimal Supersymmetric Standard Model p| looks nowadays as a very promising 
candidate for the physics beyond the SM. It can be probed experimentally through the high 
energy production of superpartners. Since no supersymmetric particles have been found 
so far, one can limit the parameter space of MSSM using the supersymmetric radiative 
corrections to the low-energy observables. The limits from the precise measurements of 
EDMs of the neutron, electron and heavy atoms are known to be a very good source of 
information in this respect [^. Namely, the one- loop contribution to the electric dipole 
moments of quarks allows one to put stringent limits on the CP- violating phases in the soft- 
breaking sector, < few x 10^'^ |0]. Alternatively, the squark masses have to lie in the TeV 
range if we believe that the phases are not suppressed. This dilemma is often referred to as 
the Supersymmetric CP problem. 

To avoid this problem from the very beginning, one can postulate that the soft-breaking 
sector as well as the /x-term in the superpotential are CP-conserving. This choice of parame- 
ters is often suggested by supergravity @]. In addition to that, it is reasonable to assume the 
universality in the mixing of left- and right-handed scalar quarks to avoid unwanted FCNC 
at low energies. As a result, the only source of CP- violation resides in the Yukawa couplings 
and can be described at the electroweak scale by the Kobayashi-Maskawa (KM) phase. In 
principle, the prediction for the neutron EDM in MSSM can be different from the SM one 
due to the specific supesymmetric contributions and the KM-related phase which appears in 
the squark mass matrix |Q. 

The order-of-magnitude estimates of EDMs [^ 0], handling the main dependence on 
Yukawa couplings and mixing angles, suggest that the EDM of the neutron does not exceed 
10^'^^e ■ cm. This means that in this type of models the contribution from supersymmet- 
ric loops never exceeds the long-distance contribution to EDM from the usual SM jO. To 



improve the order-of-inagnitude estimates, one has to solve the nonhnear renormahzation 
group equations which is difficult to do analytically. Surprizingly enough, the attempts 
to solve these equations numerically by plugging in the values for quark masses and mix- 



ing angles measured at the electroweak scale give somewhat bigger results p, |^, |T0[. The 
discrepancy between [j^ and 0, for example, constitutes six orders of magnitude so that 
it definitely requires an explanation. To resolve the apparent discrepancies and find the 
connection between the two approaches, we solve the renormahzation group equations an- 
alytically, assuming that the renormahzation group coefficient, t = (47r)~^ log(A^/M^), is 
small. In this case the equations are solvable by iterations and an analytical result for EDM 
can be obtained. It turns out that a finite result arises already after the ffist iteration, so 
that (In is proportional to the ffist power of t. This result may be viewed as the intermediate 
step between parametrical estimates and numerical calculations. The fact that t is not small 
can change the final answer somehow, but not by many orders of magnitude if we believe 
that all Yukawa couplings remain in the perturbative regime in the whole interval from Mw 
to A. In any case, this question is more of methodological interest because of the expected 
smallness of the result. 



As it was shown in [^, |T2|, the existence of a possible fourth generation in the framework 
of the usual SM gives a significant enhancement to EDM. The result, however, is still too 
small to be probed experimentally. With the natural assumptions about mixing angles it 
does not exceed 10~^^e ■ cm |T^. In this paper we consider the MSSM with four generations 
and show that the value of the EDMs is significantly enhanced. This may give nontrivial 
constraints on the parameters of the model, complementary to those coming from the direct 
search of particles from the fourth generation and from the analysis of the renormahzation 



group evolution of Yukawa couplings |jT3|, |l4| 



We organize the paper as follows. In the next section we obtain new, apparently stronger, 
limits on CP-violating phases in the soft breaking sector mediated by the chromoelectric 
dipole moment (CEDM) of the strange quark and present the analytical calculation of EDM 
in the constrained MSSM with three generations. In Section 3 we estimate the EDMs in the 
MSSM with four generations. Our conclusions are sumarized in Section 4. 



2 EDM in MSSM with three generations 

We start by writing down the structure of the CP-odd effective Lagrangian dini< 6 at the 
scale of 1 GeV which gives a major contribution to the EDM of the neutron: 

(1) 
where the sum over i runs over hght flavours, m, d and s. All other possible CP-odd operators 

dim=6 are irrelevant (see, for example, Ref. |11[). In SM we have to take into account also 

the so called long-distance contributions: the combination of two flavour-changing operators 

saturated by all possible hadronic states 0. 

Here we concentrate ourselves on the calculation of EDMs, (ij, and colour EDMs di. The 

relative meaning of these operators for the electric dipole moment of the neutron is quite a 

controversial subject. To get the final result, one has to know how to calculate the EDM of 

the neutron in terms of di and d^. This is a very difficult task if one takes into account the 

strong dynamics at large distances. To simplify this task, the naive nonrelativistic quark 

model is used very often. It relates rfjv with d^ and d^ in a simple manner: 

4 1 

dN = -dd - -du (2) 

At the same time, the chromoelectric dipole moments as well as the electric dipole moment 
of s-quark are believed to be suppressed. The importance of the information provided by 
the EDM data suggests, however, that the naive formula based on the non-relativistic quark 
model must be abandoned and replaced with more elaborated approaches handling long- 
distance dynamics. 

Let us suppose for a moment that we have the approximate universality in the sector 
of EDMs and CEDMs, for example, dd/{emd) — ds/{ems) ~ a, dd/{fnd) — ds/{ms) — b 
and a ~ b/gg. In a recent work [|l5l, the contribution of different quark EDMs to the EDM 



of the neutron were treated using the proton spin experimental data implying a non-zero 



content of the strange quark spin in the nucleon. In other words, the authors of Ref. ||15 | 
assume the identity between axial, (Ag)^, and tensor, {Aq)T, charges. (The definitions for 
these quantities are: N-f^-f5N{Aq)A = {N\q-f^-f5q\N); Na^,N{Aq)T = {N\qa^,q\N)^.) As a 
"* The EDM Lorentz structure is reducible to Ot using the identity 20-^1,75 — ief^^^apf^afi 



result, a sort of cancellation between strange and down quark contributions to neutron EDM 
was observed. This leads to apparently milder limits on </> |jl5|. It is clear, however, that 
the axial and tensor charges correspond to different structure functions and do not have to 
coincide [0 . It is especially true for the tensor charge of the s-quark in the neutron. Repre- 
senting basically the sea quark contribution, the matrix element of the strange quark should 
be especially sensitive to the property of the operator with respect to charge conjugation. 
Since C(g7^75g) = +57^75? and C{qa^yq) = —qcTfiul, we deduce that {As)t < (As)yi. The 
calculation in the instanton-inspired model shows, in fact, that {As)t ^ (As)^ [|l^. Taking 
into account the results of the tensor charge calculation on the lattice [l^, within QCD sum 



rules [jT9[ and in instanton-inspired model |T^ , we conclude that in the quark EDM channel 
the naive formula (|^) is roughly correct and, at the same time, the electric dipole moment 
of the s-quark is unimportant providing no conspiracy among different flavours. 

In our opinion, the question of conspiracy among the EDMs of different flavours is not 
relevant in MSSM on account of large contributions coming from the chromoelectric dipole 



moments of quarks. We would like to apply here the results obtained in Refs. ||2y, 0, ^ ^ 
in the context of chiral perturbation theory and QCD sum rules. If the above mentioned 
universality is held, the contributions from the chromoelectric dipole moments to the neutron 
EDM tend to dominate over the quark EDM contributions. The resulting estimate is given 
by 0,123,^: 

dN ^e{0.7dd + 0.1ds) (3) 

and normally edi > di due to the large Qs factor in di. Moreover, the second term in 
the estimate (|^) tends to dominate due to the large irts. Simple arguments favoring the 
suppression of the strange quark contribution in the quark EDM channel do not work for 
CEDM contributions. (For a recent discussion on the strange quark matrix elements in the 



nucleon see, for example, Ref. [24 



It was shown recently that this operator produces a large neutron EDM in the su- 
persymmetric 5*0(10) model I^Sj. In the Minimal Supersymmetric Standard Model with 
CP-violating phases in the soft-breaking sector, the gluino exchange diagram leads to the 
following CEDM of a quark: 

di^QsT;^ -^ F{m\ Jm^) (Asin^^ + /itan/^sin^s), (4) 

where (^a = hiv{Am\,J and (f)B = Ihi(/^^a3) ^^^ the specific CP-violating phases which 



appear in this diagram. The invariant function in Eq. (^) contains the dependence of the 
super symmetric masses in the loop. It is normahzed in such a way that -F(l) = 1. Eq. (D 
gives us the hmits on the CP-violating phases considerably stronger than those coming from 
the EDMs of quarks. Taking m^(100 GeV) = 80 MeF, /x ~ A ~ m and using (D, (|) and 
experimental limits we obtain the following constraint on the CP-violating phase: 

^„„(I55^)%8.10- (5) 

When obtaining the limit (^) we have assumed that both phases can be constrained inde- 
pendently, i.e. no fine tunning among A, /itan/3, (pA and (pB- It is interesting to note also 
that the function F reaches its maximum, F ~ 2.7 for m\^ ^ 0.2 m, so that the constraints 
on the phases in this case are almost three times stronger. 

These limits are very restrictive and from naturalness reasons we assume further on that 
the soft breaking sector does respect CP-symmetry. In what follows, we calculate SUSY con- 
tributions to the coefficients in (|I|) due to the presence of CP-violation in Yukawa couplings. 

We take the superpotential of MSSM in the following standard form: 

W = UY^QHu + DYaQH^ + EY^LH^ + piH^Ha (6) 

The soft breaking sector at the unification scale comprises flavour-blind scalar mass terms: 
m\J\Hd\^ + m]jJ\Hu\^ + mlU^U + mlb^b + nilQ^Q + H.c; (7) 

so called A-terms proportional to Yukawa couphngs 

A(uiY^QLHu + biYuQLHd + b^RY^QLHd) + H.c; (8) 

B term, 



B^H^-Hd + H.c- (9) 



and the gaugino masses 



E ^X^iR^^L +H.C. (10) 

1=1,2,3 

In formal language, the absence of additional CP-violating phases means the following: 

A = A\ B = B*, fi = /i*, m,, = m* (11) 



The squark mass matrices at the unification scale 

m^ + MtM„ (A + /icot/5)Ml 

Mi=\ ^ I (12) 

M« ( A + /i cot /3) ml + M^Mt 

m^ + MJ^Md (A + /i tan /5)M)j 
M|=| I (13) 

Mrf(A + /i tan /5) m% + M^M^ 

do not lead to any flavour- changing effects, not mentioning CP-violation, since they are 
diagonalizable in the generation space by the same bi-unitary transformation as quark mass 
matrices. As a result, quark-squark-neutralino or quark- squark- char gino interactions do 
not develop flavour- changing vertices. The terms bilinear in M„ and M^ in the squark 
mass matrices (|T2|) and (p!3| ) originate from the F-term in the superpotential. Since they 
are not related to the renormalization group running from large scale down to M^/, we 
shall refer to them as the tree-level terms. The renormalization group evaluation of the 
soft-breaking parameters and Yukawa couplings down to the electroweak scale indroduces 
flavour-changing entries and brings, for example, the dependence of Y^ in M? [^. At the 
scale of the electroweak symmetry breaking the squark mass matrices can be parametrized 
as follows: 

M? = M^l = (14) 

where different blocks have the following meaning: 

MIll = ml, + MiM^ + c,YiY^ + f^YlY, + C2{YlYuf + f2{YlYdf + 

+4{YtY„,YiY4 + ... 
M'^LR = iM'^nL)^ = iA,+ficotf3)Mi + AihY\Y, + hYiY^ + ...)Mi 
Minn = ml, + MMu + hY^Yl + h{Y^Yif + l'^Y^Y\Y,Yi + ... (15) 

and similarly for the blocks of M'j^. The explicit form of all coefficients can be found in 
the Appendix. In the expression ([15D all quark masses and Yukawa couplings are taken at 
the electroweak scale. The dots stand for other terms in the series of the increasing power 
of Yukawa couplings and renormalization group factor t. The subscripts of the coefficients 
Ci, C2, /i, ... denote the smallest power of t in which these entries to the mass matrix can 
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arise, i.e. /i ~ t^, /2 ~ t^, etc. Please, note also that C2, /2, etc. are not the contributions 
from the two-loop beta function. Finally, Ai, toqi and m^j^ are soft-breaking parameters 
renormalized by gauge interactions and by Yukawa interactions effectively conserving flavour, 
i.e. proportional to Tr[YuY^]. 

As a result, one-loop diagrams with chargino and gluino inside the loop. Fig. 1 and Fig. 2, 
can develop imaginary parts and thus lead to a nonvanishing EDM. This effect, however, 
is severely suppressed by the combination of Yukawa couplings and mixing angles. The 



extraction of the imaginary part gives almost the structure of the Jarlskog invariant pSf and 
leads to the following simple estimate Q: 

dd ~ e(loop factors) x \lm^[yidYtYlYX^^ e(loop f.) x Jcp'^Y^^Y^YI (16) 

where Jcp = Im(l^^V^^l/'^V^^); rn? is the characteristic momentum in the loop of the order 
of supersymmetric masses; "loop factors" denotes numerical coefficients reflecting the loop 
origin of the effect and the subscript 11 denotes the projection on the initial flavour d. Even 
with optimistic expectations for the numerical coefficients in Eq. (0), the result can hardly 
exceed 10~^^e ■ cm for tan/5 ~ 1. At the same time pure gluonic operators, as well as the 
EDM of the u-quark, are further suppressed by Yukawa couplings coming from down quark 
sector. 

To resolve apparent discrepancies with the numerical calculations, we solve the renor- 
malization group equation analytically by iterations retaining the smallest power of t. This 
procedure is very simple; it brings analytical answers for all operators in Eq. (^ and at the 
same time it gives the possibility to check all numerical calculation in an easy way. 

In both papers quoting the estimate (|T6|) , Ref.0 and Ref. 0, the effect is claimed to arize 
in the third or fourth order in renormalization group coefficient t. Surprizingly enough, we 
arrived to somewhat different conclusion. The chargino exchange diagram gives a nonvan- 
ishing dd in the lowest possible (first) order in t. Clearly, this contribution is related with 
the Mu-dependence of the u-squark mass matrix (|1^ coming from the tree-level and is not 
associated with t. The gluino exchange diagram yields EDMs in the down-quark sector in 
the third order in t. 

Chargino contribution 

The chargino exchange diagram which contributes to EDMs and CEDMs in the down-quark 



sector contains the u-squark line. The result of the renormalization group running, up to 
the first power in t is summarized by the set of formulae displayed in the Appendix. Here 
we give the truncated form and list only relevant entries in ([T5|) with F^^-dependence: 

/ ml. + Mt.M,, + /lY^Yrf {A, + fi cot (3)Ml + Ah^YlYMl 



Mi 



ml, + M^M, 



1^ (Ai + /i cot /?)M„ + AhiM^YlYd 



m, 



ui 



M„Mt 



(17) 



where 



/i 



log(A^/m^ 
167r2 



mQ + mu + mjj^ + A' 



h 



log(A^/m2) 



Assuming the equality between m^u and m,Q at the high-energy scale, m,u = uiq = m, we 
can conclude that the difference of m,ui and rrigi is also proportional to t. 

The next step could be the expansion of the squark propagator in AlSAu and M^ and the 
extraction of the CP- violating part from the string of Yukawa couplings and mass matrices 
0]. This procedure is completely justified for the charm quark mass since m^, Am,^ <^ m,"^ 
and seems to be ill-defined for the top flavour since we expect m,f, Arrit ~ m"^. Fortunately, 
there is no need to expand the squark propagator in terms of m^. Let us take the first term 
of expansion in A. Then the resulting flavour structure of the squark line has the following 
simple form: 

1 



Ini 



p2 - m^i 



^/lY^Y. -^^+^^^^f Mt ^ ,M. + 



1 



p2 - m^^ 



-AhYlYMi- 



1 



mIm. 



-JcpmlmtY,' 



pz 



m 



(71 



m,„mI 



M,„ 



/i (A + /i cot (3) - Ahi (m^i 



11 



Uli 



(19) 



(p2 — m?Y{j)^ — m^ — mfy 
where we took the difference of the two soft-breaking masses to be much smaller than the 
masses themselves and neglected m,"^ and m,"^ in the denominator. It is easy to see that the 
term proportional to hi cannot develop any CP-violating part unless m^m ^ mgi. However, 
since we are interested only in the first order in t, this contribution can be neglected being 
of the order t"^. In order to get the exact dependence of rrit, one has to sum up the series in 
{A + fi cot P)m.t, i.e. to take into account large i^ — iji mixing. 

This complication as well as the effect of wino-higgsino mixing are irrelevant for the ques- 
tion under study and we would like to neglect them. Assuming also for simplicity that char- 



2^ 



/) 



gino mass is equal to the scalar quark mass m and taking /i ^ — 3?Ti^(167r^)^^ log(A^/?7i 

after the trivial integration, we arrive at the following simple result: 

1 1 ,5 ,log(A^/m^) ^ mdm^mfYf^ A + II cot 3 ^ 
dd ^ -j7r^i^ed-e^) ^\ ' ' Jcp ^y^^ for m > m, 

1 1 ,3 ,log(AVm2) rridmlY;^ A + ^x cot 13 ^ .^^. 

Here we keep explicit the loop origin of Y^ and the tree-level origin of mlmj. In a more gen- 
eral case with different supersymmetric masses, chargino mixing and realistic rrit-dependence 
taken into account, the numerical coefficients and the dependence of mf and m'^ should be 
substituted by some more complicated invariant function. Its exact form is beyond the scope 
of our interest. For the chromoelectric dipole moment we have a similar result: 

~ ^3 1 log(AV"^^) J mamlmfY^ A + fi cot p 

dd ^ ^TTT^ T7^^ "^CP r for m>mt 

~ 3(73 1 log(AVm2) mdmlY^ A + /xcot/3 

^-^ - ~{KT7n T7r~2 -^cp 1 form<mi 21 

Numerically, the EDM of d-quark (pOf ) constitutes 



^-34 /tan/3 



3 



|4|~6-10-^M^1 e-cm (22) 

were we took Jcp^2x 10-^ tan/3 > 1, (levr^)-! log(AVm2) ~ 1/3 and m~A~100 GeV. 
It should be noted here that all Yukawa couplings and quark masses are taken at the scale 
of 100 GeV so that the result is further suppressed by the QCD evaluation of m'^ml from 
1 GeV to this scale. Even taking into account the enhancement associated with the CEDM 
of the strange quark we cannot obtain the EDM of the neutron bigger than 10~^^e ■ cm. 

The numerical enhancement of dd, up to the level of 10~^^e ■ cm, observed in Ref. |0 is 
nothing but an artifact resulting from a quite surprizing choice of the parameters m^scaiar = 
100 GeV and A = 5 TeV, so that A/m,scaiar = 50. (This was noted also in Ref. [l^.) We 
believe, however, that this choice of parameters is highly unnatural if not strictly forbidden. 
Here we would like to show how it could change the result numerically. In this case, the 
expansion in Am^c/m^^^^i^^ is not associated with any numerical smallness since this factor is 
now roughly of the order 1. In consequence, the result would not be suppressed in m^ at all, 
leading to an enhancement of the order 10^. 

Gluino contribution 

It is clear that the result can arise here only in the third order in t. The simple solution for 

9 



the down squark mass matrix, analogous to ( pTSD is not sufficient, though. Indeed, it can be 
shown that the flavour-diagonal projection of the combination 

^ -(A' + /icot/? + /iiYtYjM),- ^ 



d „ „ ,,>2 



P' - ml, - (^ + c,)YlY, - fMYu p' - ml, - (^ + h)Y,Yl 

(23) 

does not develop any CP-violating part at all. This shows that the simple estimate presented 

in Refs. @] and based on the quadratic anzatz of the mass matrices are incorrect since the 

exact result in this approximation is just zero in any order in t. This situation resembles the 

cancellation of the two-loop contribution to the EDM of a quark in the SM. The combination 

of different diagrams, nonzero by themselves, vanishes after the complete summation over 

all flavours ^6 . 



To be rigorous, however, we have to go further and calculate quartic combinations of 
Yukawa couplings in the squark masses and A-parameter proportional to the second power 
of t, i.e the coefficients C2, d2, etc. It can be shown also that the t'^-order entries in the quark 
mass matrices contribute to EDMs in the next (fourth) order in t and thus are irrelevant. The 
calculation is very simple and all relevant coefficients are listed in the Apendix. Omitting 
all intermediate steps, we would like to quote here the final result: 

8A{A + /itan/?) - (Sm^ + A^ " 



dd = ed-mdm^Xc^t Jcp ^ r 



-A-(A + ,t.n,3)- ^^^, 

(24) 



It is obtained with the same simplification as Eq. (|20|) and for rrit ^ m. Numerically, for 
relistic values of t and soft-breaking parameters, this result does not differ very much from 
Eq. (p2D and finally we conclude that the EDM of the neutron in constrained MSSM with 



three generations does not exceed the nonsupersymmetric SM predictions. 

3 MSSM with four generations and neutron EDM 

The extreme smallness of the result (^) kills any hopes of getting any observable effect in 
the model. Moreover, it never exceeds the non-supersymmetric SM result. The situation is 
somehow different in the constrained MSSM with four generations of fermions. This model 
acquired significant attention in recent years, since its parameter space is already severely 



restricted by the existing data |13, 14]. The detailed analysis of the renormalization group 



equations in Ref. [0 showed that the condition of perturbative evolution for all Yukawas 
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requires the masses of the fourth generation quarks to be under 200 GeV. The similarity 
between rrit and rrit' suggests that the third and fourth generations are strongly mixed, pos- 
sibly with the angle Vtb' ~ C>{^) or even ~ 0{1). This means that possible CP-odd invariant 
combination of angles involving second, third and fourth generations, lm.{V^lV^fj,V*jj,V^^) is 
bigger than Jcp in SM by a factor of the order A"^ — A^^ (See also Ref. [0). 

Clearly, this is not the main source of enhancement for CP-odd observables in the model. 
A tremendous enhancement is associated now with the change of masses and Yukawa cou- 
plings. Instead of being proportional to m'^m^Y^, all the results for EDMs contain now the 
factor ■mfmf,{mf—mf,)Y^. Obviously, this factor is comparable with m^caiar ^^^ independent 
from all "light" flavours: u, d, s, c and b. 

The detailed calculation of EDMs in the constrained MSSM with four generations is 
more difficult than in the conventional three-generation model. The main complication now 
is that not only rrit, but also mf and rUb' are comparable with m, so that we have to 
hold the exact dependence of these masses. The language of mass insertions is completely 
inadequate now. Instead, we have to perform the analysis in the mass eigenstate approach. 
In the limit of small t it is still possible to obtain the result in a closed analytical form. 
Bearing in mind, however, the big degree of uncertainty related with unknown masses and 
mixing angles, we believe that simple estimates are sufficient in this case. In particular, we 
can say that the renormalization group evolution of the Yukawa couplings and soft breaking 
parameters from the unification scale introduces at the electroweak scale CP- violating phases 
of the order Im(V^^*Vjft,Vc^,VcJY^^Y(?(Fj^ - Y^^)Y^. Thus, we can use the estimates of the 
EDM in the presence of CP- violating phases with the simple substitution: sin{(j)A{B)) ^ 
lm.(y-i*gV-i-fj,V*f,,V^g) . Then for the chromoelectric dipole moment of the s-quark we have: 

ds - 9s^ ^ MV;sVt,V:,Vj (25) 

The resulting estimate for the electric dipole moment of the neutron now is: 

dM ~ MV:yt,'V:,Vj (^^^)' 2.5 ■ 10-^=^6 ■ cm (26) 

where we took the supersymmetric mass to be of the order of 200 GeV. The overall enhance- 
ment factor in the MSSM with four generations in comparison with the three-generation 
case is enormous, being roughly of the order 10^. As a result, the EDM of the neutron is 
predicted in this model to be at the measurable level since we expect lin(yt*Vtfj,V*fj,V^g) ~ 

11 



A^ — A^ ~ 10 ^ — 10 ^. The current experimental limit on the EDM of the neutron translates 
Eq. ( p6D into the following constraint on the CP-odd combination of matrix elements: 

/200GeV\2 



lra{V,:V,,,V:,,Vji-^^^j < 4 ■ lO^^ (27) 



It is instructive to compare MSSM with four generations and nonsupersymmetric SM 
with the same number of generations. It turns out that in SM with four generations the 
result arises at three-loop level and these loops are not supported by the large logarithmic 
factors. It means that the phase space factor in the denominator is large, confining neutron 
EDM to lie within 10~^^e ■ cm fn\, O]. From the point of view of the parameters defined at 



the electroweak scale, the EDM in the MSSM with four generations arises just at one-loop 
level and therefore is considerably larger than in the non-supersymmetric case. 

4 Conclusions 



We have shown that the limit on the colour EDM of the strange quark obtained in pO| 



implies that the constraints on the phases (pA and (pB in the soft-breaking sector of MSSM 
are stronger than those from the quark EDM channel 0]. This strengthens the so called 
supersymmetric CP-problem and may have an impact on some applications of the specific 
SUSY CP- violation which use the phases on the edge of the previous limits [0. The 
solution to the supersymmetric CP-problem may be related, see for example Ref. |^, with 



the dynamical suppression of gaugino masses and/or A parameter. Alternatively, all phases 
can be chosen zero at the unification scale 0]. 

In the MSSM with the CP-violation coming from the Yukawa matrices simple order-of- 
magnitude estimates and the numerical calculations are known to disagree by many orders of 
magnitude. We demonstrate that analytical results are possible for EDMs in the limit of small 
t. The chargino-exchange diagram yields non-vanishing EDMs for the down-quark sector in 
the lowest possible (first) order in t. The result shows explicitly the suppression factor 
associated with Yukawa couplings and basically confirms the order-of-magnitude estimates 
10, 0. It provides the possibility to check all numerical calculations by lowering the high- 
energy scale. We believe, however, that even for t ~ 1/2, corresponding to A ~ 10^^ GeV the 
results cannot differ from ours by many orders of magnitude. The contrary would be possible 
only if a kind of singularity exists in the solutions of the renormalization group equations at 
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some point between My/ and A. We reject this possibility and therefore conclude that the 
EDM of the neutron in the constrained MSSM does not exceed its SM prediction 0. 

The four-generation modification of the MSSM can be viewed as an interesting gener- 
alization of the conventional MSSM and at the same time it can be checked and/or ruled 
out experimentally in the nearest future |T^, |T^ . We believe that the similarity between rrit 
and irtt' suggests also non-zero, presumably large, mixing between the third and the fourth 
generations. At the same time, the relevant combination of masses is much bigger in this 
case and the enhancement factor is given by (m^,m^)/(m^m^) ~ 10®. This leads to a large 
EDM of the neutron at the measurable level of 10~^^e • cm. 
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Appendix 

A) Assuming that the renormahzation group coefficient t is small, the renormalization 
group equations can be solved analyticaly by iterations. In this case, the coefficients which 
appear in equation (^5|) are found to have the following explicit form: 

h^ = -t + 0{f) + ... (28) 

k^ = -3t + 0{t^) + ... 

These quantities arise after the ffist iteration and are the relevant entries to the u-squark 
mass matrix for the calculation of chargino-exchange contribution to quark EDM. 

Since gluino-induced quark EDMs appear only in the third power in the renormalization 
group coefficient t, we have to perform also the second iteration. The results for C2, /2, c'2, 
I2 and ^2 are: 

C2 = 3t^A^ + 0{t^) + ... 

/2 = 3t^A^ + 0{t^) + ... 

4 = t^A^ + 0{t^) + ... (29) 

I2 = 6t^A^ + 0{t^) + ... 

1'2 = e[ml^^+mX+ml^^+3A^]+0{t') + ... 

These expressions are obtained starting with diagonal soft-breaking mass parameters at the 
unification scale; the index refers to the parameters at that scale. In order to obtain the 
coefficients for the entries in the d-squark mass matrix, one has only to interchange the 
indices u and d. 

B) We would also like to present the complete solutions up to the second order in t for 
the renormalization group equations associated with the soft-breaking parameters. We used 
the notations and the one-loop /3-functions given in [^ . 

+e A^[?,{YlY^.,f + 3(Yt^Y,j2 + {Yt Y„„ Yt^Y.J] 
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m 



U2 



m^ 



ml + a"Y„Yt + 
mi + a'^Y.Yi + 



'd2 



+]f[l2AlY,,YlY,,Yl + {AAl + 6^)Y,,Yt^Y.,Yi 



d2\ 



•■U2 ^ U2 



Ida — Ida 



A^ + a^"Yt Y„ + fo'^" Y;^Y, 
Aa + a'^'YlY, + 6"^Yt Y„ 



In the above formulae we used the following notations: 







(30) 



6" 



-t m. 



% -r "^//^ 



"^i+"^7f. +^1 + T (mir,^ -mir)dt + 






+^'[Pqo + P^o + ^^'^o] - ^I2m^„ + 2m4, + 2m' + 3 A 



Z 

«,d^ _2^r^2 , 2 ^ 2 _^ ^21^ 1/2 _^2 Nj^^ 



+ 7:t' 



Jo + ^^«,do - 2/?«,do["^Qo + "^l^^o + "'L.dn + ^ 



6"''^ = 2 fm^ + mi, + m?. + A^l 



6'^" 



(31) 



Finally, A^^, mgj, m„2 ^^^ m^j are soft-breaking parameters renormalized by gauge in- 
teractions and flavour-conserving Yukawa interactions. The indices 0, 1 and 2 refer to the 
iteration to be considered for the corresponding quantity. The functions r^, p and /? are 
scalar quantities involving gauge couplings, gaugino masses and Tr[Yj- Yj]. Their exact form 



can be extracted from Ref. ||29||. Using these solutions, one can identify all the coefficients in 
equation ([T5|) up to the second power in t. 
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Figure captions 

Fig. 1. Quark self-energy involving chargino exchange, generating EDM (CEDM) of quark 
in the external electromagnetic (color) field. 

Fig. 2. Quark self-energy involving gluino exchange, generating EDM (CEDM) of quark in 
the external electromagnetic (color) field. 
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